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^ Abstract 

We develop Bogoliubov theory of spin-1 and spin-2 Bose-Einstein condensates (BECs) in the 
presence of a quadratic Zeeman effect, and derive the Lee- Huang- Yang (LHY) corrections to the 
ground-state energy, pressure, sound velocity, and quantum depletion. We investigate all the phases 
^ of spin-1 and spin-2 BECs that can be realized experimentally We also examine the stability of 
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2 each phase against quantum fluctuations and the quadratic Zeeman effect. Furthermore, we discuss 



a relationship between the number of symmetry generators that are spontaneously broken and that 
of Nambu-Goldstone (NG) modes. It is found that in the spin-2 nematic phase there are special 



s 

Bogoliubov modes that have gapless linear dispersion relations but do not belong to the NG modes 
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I. INTRODUCTION 



The Bogoliubov theory of weakly-interacting Bose-Einstein condensates (BECs) [T] has 
served as an indispensable tool in diverse subfields of physics. For a scalar BEC of bosons 
with mass M, particle- number density n, and 5- wave scattering length a, the ground-state 
energy (GSE) of the system with volume V is given by 



where the first term on the right-hand side is the mean-field energy, the second term gives 
a nonperturbative correction to it which was first derived by Lee, Huang, and Yang (LHY) 
[21 |3j, and the higher-order terms were discussed in Refs. |4]-[6j. In the present paper, we 
discuss Bogoliubov theory and LHY corrections of BECs with spin degrees of freedom in 
the presence of a quadratic Zeeman effect. 

The Bogoliubov theory of spinor BECs has been discussed extensively over the past 
decade. The spin-1 Bogoliubov spectra have been derived in Refs. |7l-[T0] up to the linear 
Zeeman effect. In Ref. [llj, the same problem is discussed from a field-theoretic point of 
view. The effect of the quadratic Zeeman energy on the spin-1 BEC has been discussed in 
Refs. [121 US]- The spin-2 BEC has been examined in the absence of an external magnetic 
field in Ref. [M] and up to the linear Zeeman effect in Ref. [15j. However, little attention 
has been paid to the GSEs. 

In this paper, we develop a systematic renormalization procedure and derive GSEs, pres- 
sure, sound velocity, and quantum depletion up to the LHY corrections. The LHY correc- 
tions have been measured for a scalar BEC p!6l \T7 \ and for a two-component Fermi gas [TS] 
by using methods to enhance quantum fiuctuations. The experiment in Ref. [16j utilized 
a strongly correlated system in an optical lattice, while the experiments in Refs. pTl [T8] 
amplified the coupling constant by means of a Feshbach resonance. Our analysis takes into 
account the quadratic Zeeman effect that is of great importance under many experimental 
situations in which the linear Zeeman effect can be ignored. Because the sign of the quadratic 
Zeeman term q can be manipulated experimentally [19j, both cases of positive and negative 
q are analyzed. It is shown that except for the ferromagnetic phase the LHY correction in 
a spinor BEC is affected by the quadratic Zeeman effect and that it can be measured by 
making strongly correlated systems or by controlling the external magnetic field. 
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The order parameter of a spin-2 nematic BEC in the absence of an external magnetic 
field depends on an additional parameter, ry, that is not related to the symmetry of the 
Hamiltonian but describes the degeneracy between the uniaxial and biaxial nematic phases. 
As pointed out in Refs. [20l [2T], however, quantum fluctuations induce a quantum phase 
transition between the two phases, lifting the degeneracy. We show that the quadratic 
Zeeman eflFect with q < causes the dynamical instability in the uniaxial nematic phase, 
whereas it leaves the biaxial nematic phase stable. That is, the uniaxial nematic phase is 
unstable against an inflnitesimal negative quadratic Zeeman effect in the thermodynamic 
limit. Conversely, the quadratic Zeeman eflFect with q > makes the biaxial nematic phase 
dynamically unstable while it leaves the uniaxial nematic phase stable. However, it is pos- 
sible to stabilize both of these phases for nonzero g in a flnite system. We will show this for 
the case of a spin-2 BEC. 

The Bogoliubov theory predicts massless modes, which can be interpreted as Nambu- 
Goldstone (NG) modes associated with spontaneous symmetry breaking. To elucidate this 
point, we discuss a relationship between the number of symmetry generators that are spon- 
taneously broken and the number of NG modes [22j. We apply the relationship to spin-1 
and spin-2 BECs, and point out that for the uniaxial and biaxial nematic phases there exist 
the Bogoliubov modes that have gapless linear dispersion relations but do not belong to the 
NG modes. 

This paper is organized as follows. Section H formulates the problem, and describes 
the low-energy Hamiltonian and Hartree-Fock approximation of a spin-/ BEC. Section I 
n discusses the problem of divergence of the GSE and how to remove the divergence by 
renormalization of the coupling constant. Sections IV and V examine the mean-fleld phase 
diagrams and Bogoliubov theory of spin-1 and spin-2 BECs, respectively, in the thermody- 
namic limit, and derive the Bogoliubov spectra and LHY corrections. Section VI discusses 
the relationship between the number of symmetry generators that are spontaneously broken 
and the number of NG modes in spinor BECs. Section VH provides the summary and con- 
cluding remarks. The detailed derivations of the GSEs are described in Appendix A, and 
the properties and equation numbers of the physical quantities in each phase and notations 
are listed in Appendix B. 
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II. FORMULATION OF THE PROBLEM 



We consider a system of spin-/ identical bosons with mass M that undergo an 5-wave 
scattering subject to periodic boundary conditions. As in most experiments done in spinor 
BECs, we consider the case in which the hnear Zeeman effect can be ignored. Let ^^rn(x) 
(m = — /, — / + !, ...,/) be the field operator of a boson at position x with magnetic quantum 
number m, where we assume that an external magnetic field B is applied in the z direction. 
Then, the low-energy eflFective Hamiltonian of a spin-/ BEC is given by 

H = Hke + Hqz + V, (2) 

where 

^KE = |rfx ^L(X)( - (3) 



is the kinetic energy, 



Hqz = qm^ J d^K ^^|^(x)^^^(x) (4) 



is the quadratic Zeeman term, and 

2/ F 

^ = E f E {fmfm'\FM){FM\ff,ff,') j dx ^'L(x)^^L(x)^^^(x) V(x) (5) 

F=0 M=-F 

is the interaction energy. The strength of the quadratic Zeeman term is given by q = 
{qIIbBY / Ei^i^ where g is the Lande ^'-factor, jiB is the Bohr magneton, and E'hf is the 
hyperfine energy splitting. In Eq. ([5]), ^i? is a bare coupling constant in the total spin F 
channel and {fmfn\FM) is the Clebsch-Gordan coefficient. Here and henceforth, repeated 
indices such as m^m'^jji^ and ji' are assumed to be summed over f^f — 1,...,—/ unless 
otherwise stated. Bose symmetry requires that the total spin F in the 5- wave channel 
is even. In fact, it follows from the canonical commutation relations of bosons and the 
properties of the Clebsch-Gordan coefficients that the terms in Eq. ^ with odd F vanish 
identically. 

We expand the field operator as 

^^^(x) = ^ V ak,.ne'*"'', (6) 

VV ^ 
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where V is the volume of the system and ak,m is the annihilation operator of a spin-/ boson 
with wave number k and magnetic quantum number m. Substituting Eq. ^ in Eq. ([2]), 
we obtain 

2/ _ F 

k F=0 M=-Fk,p,q 

where 6k h^\s?/2M. 

In the mean-field or Hartree-Fock approximation, all bosons are assumed to occupy a 
single mode, which we label as 0: 




where variational parameters are assumed to satisfy the normalization condition 
X^mlCmP = 1 and they are determined so as to minimize the expectation value of the 
Hamiltonian. By using the trial state ([s]), it is possible to classify the mean-field ground- 
state phase in a number-conserving manner p!Ol [T5l [23] . 



III. RENORMALIZATION OF THE GROUND-STATE ENERGY 

If we take into account the eflFects of quantum fiuctuations, physical quantities such as 
GSE show ultraviolet divergence. This divergence stems from use of the contact interaction, 
which gives correct results if we only consider the region |x| ^ rg, where Tq is the range of 
the interaction. This implies that the effective Hamiltonian ^ is only valid below a certain 
cutoff momentum. By introducing the cutoff, the GSE no longer diverges but depends 
explicitly on the cutoff. The renormalization of the coupling constant eliminates the cutoff 
in favor of an observable, that is, an 5-wave scattering length in the present problem. 

To examine this problem of the coupling constant in detail, let us first consider the low- 
energy scattering between spin-/ identical bosons. The scattering rate for each scattering 
channel is determined by the T-matrix T, which is the solution to the following equation 



f = V + V — f, (9) 

E — 2i7KE 

where E is the total energy and !/(£' — 2/fKE) is the two-particle Green's function. We deal 
with the quadratic Zeeman term perturbatively by assuming that it is at most of the order 



of y - the condition well satisfied in current experiments. This assumption is implicitly 
made in classifying the mean-field ground-state phases in Refs. [121 EH 128] . 

The T-matrix is related to the 5- wave scattering lengths and renormalized coupling 
constants SlS follows: 



9f 



M 



lim(k^F|T|k,F), 



(10) 



where k and k' are the incoming and outgoing wave vectors, respectively, and A; = |k| = 
|k^| because of energy conservation. In the Hartree-Fock approximation, the T-matrix is 
approximated hy T ^ and therefore 



9F = ^^=9F. 



In the Bogoliubov theory, however, this approximation does not remove the cutoff depen- 
dence because the divergence occurs at the level of second order in g^. Thus, we must 
approximate the T-matrix up to second order in the couphng constants: 

1 



T^V + V- 



-V. 



(12) 

E — 2i7KE 

The corresponding relation between the bare and renormalized coupling constants is given 
by 



k 



(13) 



where the prime on the sum means to omit the term k = 0. The diagrammatic representation 
of Eq. (13) is shown in Fig. [!} 





FIG. 1: Diagrammatic representation of the coupling-constant renormalization in Eq. (13), where 



each solid line describes the one-particle Green's function of the Hamiltonian ^/^ke and each vertex 
represents V . The first diagram on the right-hand side represents the Hartree-Fock term, and the 
second diagram represents the second-order term. 
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Next, we illustrate how this renormalization of the coupling constant eliminates the cut- 
off dependence. As we will show in subsequent sections, the GSE in the absence of an 
external magnetic field is expressed as 

E, = + E^^. (14) 



Here is the mean-field energy 



Er = (15) 



where n = N/V, and C* is a linear combination of the bare couphng constants. The term E^^ 
in Eq. ( [I4) ) describes the contribution from quantum fluctuations around the Hartree-Fock 
mean field and takes the following form: 



fc2 / 

^0 --Ji^l^l^ 



4M 

3 k 



,2 2MnC7. , 1^ AMnCi 



(16) 



where the sum over j is taken over the Bogoliubov modes, each of which describes fluc- 
tuations such as density and spin fluctuations, and Cj is a linear combination of the bare 
coupling constants. At the limit of A; ^ cx), the integrand behaves as follows: 

2 



where we substitute a renormalized coupling constant Cj for Cj on the right-hand side, which 
is correct up to the second order in the coupling constants. On the other hand, from the 



T-matrix calculation of Eq. (13) (note that gjr and correspond to C and C, respectively). 



the mean-fleld energy is calculated to give 

,2/^2 fc2 / /OA/r^/^.\2 



^° =^ + m^^ [-iPtT ■ ^^^^ 



Here the second term on the right-hand side cancels the ultraviolet part of Eq. (16). Thus, 
the cutoff dependence of E^^ is removed. Therefore, the GSE is given by 

2 ' 



^0 - -^-Jul^l^ 



2 4M 

i k 



^2 ^ 2MnCj _ / , , 4MnCj 1 f2MnCj 
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where Xj = 27r/i/(y^/^y^2MriCj), and we take the thermodynamic hmit Xj = to obtain the 
last expression. Even if we incorporate the quadratic Zeeman effect, the above canceUation 
mechanism holds as shown in Appendix A. 



IV. SPIN-1 BEC 



For a spin-1 BEC, the total spin F of two interacting bosons must be or 2, and therefore, 
Eq. ^ reduces to 

k k,p,q \ 

(1)^ ^ .-j- .-j- . 



(20) 



where = (^o + 2^2)/3, c^i^ = {92 - c/o)/3, and f^^/ = (/^^,, /^^,, /^^,) represents a set 
of the spin-1 matrices given by 



r 



V2 



/o 1 o\ 

1 1 
^OlOy 



V2 



0-10 

1 -1 
1 



, r- 



1 
000 
00-1 



(21) 



The possible phases and phase boundaries of spin-1 BEC are shown in Fig. [2] and de- 
scribed as follows |12l [27] : 



ferromagnetic phase: < 0, c^^'' < 0, and C,^ = (1,0,0), 
polar phase: (? > 0, q> —2c^^^n, and C'^ = (0,1,0), 



g<0, ^^>0, and (^' = -^(1,0,1), 



1 



(22) 
(23) 
(24) 



(i)/n o.(i)^ /^sin^ . sin^ 



broken- axisymmetry phase: Q'>0, c\ q < —2c\ n, and ( 



\V2 



, COS^, 



(25) 



where = ((?o + 2g2)/3, = ((?2 - (?o)/3, sin^ = y^l/2 + g/(4n4')). In Fig. |2| the 
shape and color of the wave function in each phase represent the symmetry of the order 
parameter. For example, the spinor of the polar phase, has a rotational symmetry about 
the z axis, since the shape and color are symmetric about the same axis. The ferromagnetic 
phase apparently does not have an axisymmetry because the color changes around the z axis. 




F P' 



FIG. 2: (Color online) Phase diagram of spin-1 BECs, where F, P (P'), and BA stand for the 
ferromagnetic, polar, and broken-axisymmetry phases, respectively. The spinor order parameters 
(C15C05C-1) of P P' are given by (0,1,0) and (1/v^, 0, l/\/2), respectively, which are trans- 
formed each other by a 7r/2 rotation about the x axis. The degeneracy of these states are lifted 
by the quadratic Zeeman effect. The thick lines represent the phase boundaries. The shape of 
the wave function in each phase is calculated in terms of polar angle and azimuthal angle 
as |^^(6>,0)p = |EmCmY7'(6>,0)p, where yr(6>,0) (m = 1,0,-1) are the spherical harmonics of 
rank 1 and color represents the phase of ^^(0,0) (see the right color gauge). The axes of the spin 
coordinates are shown at the bottom right of the phase diagram. 

However, the U{1) gauge transformation can make up for the variation of the color so that 
due to the spin-gauge symmetry, the ferromagnetic phase maintains the U{1) symmetry. 
The broken-axisymmetry phase does not possess any continuous symmetry because any 
U{1) X 5^0(2) transformations cannot compensate for the variation of the color shown. 

The ferromagnetic phase has a longitudinal magnetization, whereas the broken- 
axisymmetry phase has a transverse one that depends on q. Moreover, the broken- 
axisymmetry phase has a finite spin-singlet pair amplitude that also depends on q. In 
contrast, both spinors of P and P' have no magnetization but a finite spin-singlet pair 
amplitude that is independent of q. These properties in each phase are summarized in 
Appendix B. Since c^^^ > 0, the spin-1 ^^Na condensate is in the polar phase. On the 
other hand, since c^^^ < 0, the spin-1 ^''Rb condensates can be in any of the ferromagnetic, 
broken-axisymmetry, and polar phases, depending on the sign and magnitude of q. 



In the Bogoliubov theory, we replace operators ao,m by c-numbers V^iVoCm and keep 
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and ak/o,m up to the second order in the Hamiltonian, where A^o is the number of condensate 
bosons, which, together with nk,m = Skm^k,m5 satisfies 

No + J2J2<"^ = (26) 

k m 

where is the total number of bosons. The Bogohubov Hamiltonian of a spin-1 BEC is 
thus given bv fT0lfT2] 



(4ncf)(f)^)+giv((rr) + j; 



(1) 
L 

(1) 



nc, 



(1) 



(2Fj, • Fk + Fk • F-k + Fj, • F^_k) 



(27) 



where 



m 



(28) 
(29) 
(30) 



Here, and Fk denote the density and spin fiuctuation operators of the condensate, 
respectively. In Eq. (27), we substitute c^"^^ for c^"^^ {i = 0, 1) in the sum over the momentum 



because the Bogoliubov approximation is correct up to the second order in the coupling 
constants. 

A. Ferromagnetic phase 



For the ferromagnetic phase (22), Eq. (27) reduces to 



^eff 



+ qN + Y^ 



(6k - oSk,o + (^k - 2nc^^^^ a[ _iak -1 



+ (^^k + n[CQ + )j %,i<^k,i H ^ l%,i^-k,i + <^k,i<^-k,ij 



(31) 
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Here, the m = and —1 modes are already diagonal, and the m = 1 mode can be diagonal- 
ized by the standard Bogoliubov transformation [71 18] : 



2£'k,i 



Sk,i + 



'ek + n(c«+c«)-^k,i 



2£'k,i 



a' 



k,i' 



where E\^ i is the Bogoliubov spectrum given by 



The diagonalized Hamiltonian is 



ek 



ek + 2n(4^ + 4 



(1) , 



^eff 



^k,i^ki^k,i + (ek - q')«ko«k,o + (ek - 2n4^^)aJ._iak-i 



where 



(32) 



(33) 



(34) 



(35) 



is the GSE in the ferromagnetic phase. As can be seen from Eqs. (33) and (34), the m = 1 



mode is massless, and in the absence of the external magnetic field, the m = mode is 
also massless. For the excitation energies with m = and m = — 1 to be positive, we must 
have q < and c^"^^ < 0. For the m = 1 Bogoliubov mode to be stable, Cq"^^ + c^"^^ > is 
required. This condition ensures the mechanical stability of the mean-field ground state; 
otherwise the compressibility would not be positive definite and the system would become 
unstable against collapse. These requirements are consistent with the stability criteria of 
the mean-field theory. Conversely, if we prepare a spin-polarized state with g > 0, c^^^ > 0, 
or Cq"^^ + c^"^^ < 0, the state would undergo the Landau instability for the m = and — 1 
modes with quadratic spectra and the dynamical instability for the m = 1 mode with a 
linear spectrum. The GSE per a volume V in the ferromagnetic phase is calculated to give 



El 

V 



= qn + 



M 



-a2 1 + 



128 
15" 



TT 



(36) 



where the last term is the LHY correction. The pressure is obtained from Eq. (|36|) as 



dEl 
dV 



M 




(37) 
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and the sound velocity is given by 



1 dP UTih^na2 , 



nan 



(38) 



The quantum depletion is calculated to be 

N -No 1 v^' 



8 /na| 

TT 



(39) 

k m 

We note that the LHY corrections in the ferromagnetic phase are not affected by the 
quadratic Zeeman effect. 



B. Polar phase 



The polar phase has two spinor configurations (23) and (24), which are degenerate and 



connected each other by a [/(I) x S0{3) transformation for q = 0. However, for nonzero g, 
the degeneracy is lifted, and they should be considered as different phases. This is because 



the phase of (23) has a remaining SO (2) symmetry, whereas the phase of (24) is not invariant 
under any continuous transformation. As a consequence, it is expected that the number of 
the NG modes is different in each phase, and that the low-energy behavior is also different. 



1. Case of q > 



In the case of (23), we introduce the following density fluctuation operator ak,(i and spin 



fluctuation operators Skj^ and Skj^ around the x and y axes: 



1 . , 



In terms of them, the total Hamiltonian is expressed as 

k 



H: 



eff 



j=fx-,fy 



(1) "I 



(40) 
(41) 

(42) 



(43) 
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This effective Hamiltonian can be diagonalized by means of the foUowing Bogohubov trans- 
formations 



2E^ 



k,d 



k,d 



-k,d' 



, / ek + <? + ncS^^ + ^^kjt . , / ek + + nc^'^ - Ekjt 
Ok,j = \l 7^ -aw J + \l — -a' 



.(1) 



2^k,/t 



'^Ekjt 



k,i' 



with the result 



4ff = K + E' [^k,.sj,_/k,. + ^k,/.(S[,^^6k,/. + bifikj,)] ■ 



(44) 
(45) 



(46) 



where the Bogohubov spectra are given by 



Ek,d = \/^k(ek + 2n4^^), 
Eujt = V (ek + g)(ek + q + 2ncS'^) 



(47) 
(48) 



and the GSE by 



1. 



k 



(ek + ncj,^^ - £;k,d) + 2 (^ek + g + nc^^^ - E^^j,^ 



(49) 



The density mode (d) is massless because the U{1) gauge symmetry is spontaneously broken 
in the mean-field ground state, while the transverse magnetization {fx^ fy) modes are massive 
for nonzero since the rotational degeneracies about the x and y axes are lifted by the 
external magnetic field. 



In the limit of g ^ 0, the Bogohubov spectra (47) and (48) reduce to those obtained in 



Refs. [7|0|, and the transverse magnetization modes become massless since the symmetry 
of the Hamiltonian becomes U{1) x SO (3) and the rotational symmetry around the x and 
y axes are spontaneously broken. The conditions g > 0, Cq"^ > 0, and q > — 2c^'^^n are 
required for the Bogohubov spectra to be positive semidefinite, and they are consistent with 
the stability conditions of the mean-field ground state. 



Using the renormalization procedure of the coupling constant shown in Eq. (13), we find 
that the GSE per volume V is given by (see Appendix A for derivation) 



V 



2jl) 



n c. 



1 



16VM3 
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where ti = g/(n|c^^^|) — 1, and 
15 



dx (^x^ + t + l- ^{x^ + t){x^ + t + 2) - 



(51) 



with 01 (0) = 1. The GSE of the polar phase consists of two parts corresponding to density 
and spin fluctuation, and the latter part depends on The behavior of 01 (t) for positive 

t is plotted in FIG. |3| which shows that the GSE increases with increasing the quadratic 
Zeeman effect. 




FIG. 3: Plot of 01 (t) in Eq. (51) which describes the contribution to the ground-state energy (50) 
from spin fluctuations due to the quadratic Zeeman effect, where t = ^/(^|c^^^ 



1. 



The pressure and sound velocity are calculated as 



2jl) 



n Q 







1 + 



8VM3 



21 



n(|4')|)3[30i(ti+sgn(4'))) 
2(ti + l)0;(ti+sgn(4')))l, (52) 



and 



nc, 



(1) 



c = 



M 



TV 



\/n(c«)3 + 



2VM3 / Ic 



n(|^)|)302(ti,sgn(4^O) 



(l)^ 



(53) 



where 



02(t,sgn(cS'))) = 15<^i(t + sgn(4'0) " 12(t + l)<^;(t + sgn(4'0) + 4(t + 1) + sgn(ci'0!)54) 

which is plotted in Fig. We note that for the case of c^p < 0, the quantum correction 
to the sound velocity due to spin fluctuations vanishes at t ~ 1.4. On the other hand, the 



14 



^,(?,sgn(c!")) 



40 


Cj > U 




30 






20 
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FIG. 4: Plot of Mt,sm{c?)) which gives the correction to the sound velocity (\53h due to the 
quadratic Zeeman effect. 
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FIG. 5: Plot of 03 (t) which describes the quantum depletion through Eq. (55) 



quantum depletion is calculated to be 



k 



n(c«)3 + 2Vn(|c«|)303(ti + sgn(c«)) ) , (55) 



where 



The plot of 03 (t) is shown in Fig. [5| which shows that the quantum depletion is suppressed 
by the quadratic Zeeman effect. In contrast, the quantum corrections to the GSE and sound 
velocity are enhanced by it. Since the variations of 0i(t), 02 (t, ±1), and 03 (t) in the typical 
region of t are of the order of 1 ^ 10, the magnetic-field dependence of the LHY corrections 
is significant. However, because c^^^ ^ c^"^^ in the alkali species, the magnitudes of the LHY 
corrections are small and the quantum corrections arise mainly from density fluctuations. 
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2. Case of q < 



To analyze the case of (24), we introduce the foUowing fluctuation operators: 



. _ 1 . , 



(57) 

(58) 
(59) 



where dk,d, «k,/a;i and dk,/^ describe the density fluctuation (d) and spin fluctuations {f^, fz) 
around the x and z axes, respectively. The Hamiltonian is diagonahzed by the foUowing 
BogoUubov transformations 



2£;, 



2E, 



k,d 



-k,d' 



2E, 



'ek-g + ncS^^ -Ek,/, 



2E, 



»k,/, = \/ 7^ flk,/, + \l a 



,(1) 



2E, 



2E, 



with the result 

= K + 1] (^k,ii,/k,d + ^k,/iij,&k,/. + ^k,/iij,&k,/j, 

k 

where 



2^(1) 



(ek + nc^Q^ - Eu,d^ + (ck - g + nc[^^ - £^k,/.) 

+ (ek + n4'^ -£;k,/.) 



is the GSE, and the Bogoliubov spectra are given by 



Eu4 = yek (ck + 2n4^^). 



-E^k,/, = Y (ek - g) (ek - g + 2n4^^) , 
^k,/. = ^ek (ek + 2n4')). 



(60) 
(61) 
(62) 



(63) 



(64) 



(65) 
(66) 
(67) 
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In contrast to the case of g > 0, one of the spin fluctuation modes, (67), becomes massless. 



This is because aU the continuous symmetries of the Hamiltonian are spontaneously broken 
for q < 0. On the other hand, the transverse spin mode is massive because the rotational 
symmetry about the x axis is broken by an external magnetic fleld. For the Bogoliubov 
spectra to be real, the conditions q < 0^ Cq^^ > 0, and c^^^ > must be satisfled; otherwise. 



the state in Eq. (24) would be dynamically unstable. The GSE Eq\ pressure P, sound 
velocity c, and quantum depletion (A^ — No)/N are given by 



V 



2jl) 



= qn + 



n c, 







1 + 



16VM3 



2jl) 



n(4^))3(l + 0i(t2)), (68) 



2 (1) 



1 + 



8VM3 



2jl) 



n(c;'))3(3 + 301 (^2) - 2h4>[it2)),m 



nc, 



(1) 



M 



^ I V^(cS'^)^(15 + 02(^2-1,+!)) 



, (70) 



N-No 
N 



VM3 
3^ 



yJn{c^^^)^ + ^n{c?ni + Ut2))), 



(71) 



where t2 = —q/{nc^^^). Reflecting the difference in the Bogoliubov spectra, the LHY correc- 



tions for Eq. (24) are also different from those for Eq. (23) 



C. Broken-axisymmetry phase 



The mean-fleld ground state of the broken-axisymmetry phase is parametrized by Eq. 



(25). We deflne the following fluctuation operators: 

sin^ 



72 



(ak,i + ak-i) + cos<9ak,o, 



(72) 



^k,/^ - (Sk,l - -l) . 

In addition, we introduce the following operator that is independent of a\^^d and Sk,/^: 

cos^ 



%/2 



(ak,i + ak-i) - sin^ak,o- 



(73) 



(74) 
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In terms of them, the total Hamiltonian is expressed as 

9.-(l) . -(1)^ .rO 



TjBA 
-"eff 



Nq ^ Vn\c};' +cY') ^ Nq' 



8nc\ 



+ ^ (^ak,da_k,d + <,d«-k,d j + (^«k,f?a_k,e + <,e«-k,e 



nc, 



(1) 



gsin2^ / I 



where 



k,d'^k,6' + Q-k.e'^k.d + C(.k,dC(.-k,e + Q-k,dC(._k,i 



.(1) = ^ 



(75) 



(76) 



To diagonahze the sub-Hamiltonian for the spin fluctuation mode around the z axis, we 
perform the following transformation: 



/ek + q/2 + E^,, 



«k,/, + 



'ek+i/2-£kj..t 



2-Ek,/^ ' ' V 2£;k,/^ 

where the Bogoliubov spectrum is given by 

^kj, = Vek(ek + g). 



(77) 



(78) 



On the other hand, for the density fluctuation mode and the 9 mode in Eq. (74), we consider 
the foUowing transformations: 



where the bold symbols represent the following set of the operators 
and U{k) and V{k) are 2x2 real matrices. 



(79) 



(80) 



U(k) = - 



V(k) = - 



2C,+{k)E,{k) + Ci+{k)Xi_{k) qe^Ci+ik) sin 2^ + 



2Ci-{k)Ei{k) 
1 



-2CmW) +^1-^^)^!+^^) gekCi_(A:)sin2^ + 
2Ci+(fc)£;i(fc) + '^1+ (^)^i- (^) -9ekCi+ (fc) sin 26 - 



Xi-(k) 




-{k)Ei(k) sin 26> 








^(fc)£;i(/c)sin26' 




Xi_(/c) 




-(k)Ei(k) sin 26* 







(81) 



2gekCi+(/c)£;i(/c) sin26> 



^(82) 
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with 



Eiik) 



[(n2(c?) + 3c? f - 4n2c(^)(4^) + 2c?)]ei - in'c^ici^^ + 3c?){c? 

1/2 



(A;)ek + g2e2 (ck - 2nc[^'^) sin^ 2^ 



In Eq. (85), £'k,± are the Bogoliubov spectra given by 

-E'k,± = 



(83) 

(84) 
(85) 

(86) 



Using the transformations (79), the total Hamiltonian is diagonahzed as 

/ 



fiBA 



(87) 



where the GSE is given by 



^0 ~ 2 + 2 +Q..(i) 



8nci - k 



[ (ek + <?/2 - Ek,/J 



+(2ek + nc? - nc\'> - E^^^ - E^^_) 

^ 2n(c^^ + c?)e\i in the long-wavelength limit, £'k,- is linear and massless. 
Furthermore, the fz mode is also linear and massless for nonzero q. This is because the U (1) x 
50 (2) symmetry is completely broken in the broken- axisymmetry phase. The positivity of 
the density and spin modes is ensured if g > 0, c? < 0, c? + c? > 0, and q < —2c?n. 
Even though each mode is still massless for g = 0, in the long- wavelength limit, the spectrum 
of fz changes from linear with g 7^ to quadratic with q = 0. The physical origin of this 
change is discussed in Sec. V. Using the renormalization procedure discussed in Appendix 
A, we find that the GSE per volume V in the broken- axisymmetry phase is given by 



(1) 



Since E^ 



Eq 



nq 71^(0? + c 



V 2 ■ 2 8ci 

with ts = (f' /{nc^^^ + nc^^^y. Here 



-+j4 + « + ?^M!M£!M)ll^,,y, (89) 



(1) 



157r2^3 



15 



8^ 



dx X I ^x 



yi 



it) 



it) 



1 ^ c^t 



(90) 
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\t) = y a;4 + y,x^ + y^it) ± y^{t)x^ + y^{t)x^ + yl{t), (91) 
with y, = (c« - cS^V(4'^ + 4'^ y2(t) = (24^^)7(4^^ + cS^V - V2, ysit) = (4^^ + 



34^))V(c?) + cY'f - (4^^ + 2c^nt/c\'>, and ?/4(t) = -4(4^0^(4^^ + 30^0/(4^^ + 4^0^ + 
(cq"^^ + 30^"^^)^/ (cq"^^ + c^i''). The behavior of 04 (t) is shown in Fig. [g, where the typical values 
of t for ^''Rb is 10~^. As can be seen from Fig. |6](a), the density component of the GSE 
increases slightly due to the quadratic Zeeman effect. The pressure and the sound velocity 



(1)a2/J1) 



.(1)\3 




2.x 10 - 4jk10 = 6.X10 ' E.xlO^ 



2x10 ^ 4.X10 ■ 6 x10 - &.x!0 



FIG. 6: Plots of (a) 04 (t) and 06 (^), and (b) 05 (t). Over the ranges of t shown, they are well 



approximated as 04 (t) 2^ 1 + 135t and 06 (t) 2^ 1 + 80t + SOOOOt^. (05 (t) is given by (94)) 



are 



P 2 ^ - -^^^ + (3^4(t3) - ^hUh)), (92) 



1 + 



Vn(c?)+4^))3(/>5(t3)) , (93) 



M \ 157r2n3 
where 

05(t) = 15<^4(t) - 16t0;(t) + 16tV;'(t). (94) 



The behavior of ^5(t) in Fig. |6[b) shows that the LHY correction of the sound velocity 
decreases slightly as the external magnetic field increases. 
The quantum depletion is 

k 

+ Vn(c(nc«)30e(t3)) , (95) 
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with 

6 



^ _ / dxx\VA + Vt, + V,\ + V^,), (96) 



where Vij's are the the components of V{k) in Eq. (82). The behavior of (f)Q in Fig. [6](a) 
shows that the quantum depletion increases with increasing = / [nc^^^ +nc^^^Y . However, 
the magnetic susceptibihty of the quantum corrections in the broken-axisymmetry phase is 
smaU because the variations of 04 (t), 05 (t), and 06 (t) over the typical range of t are of the 
order of 10~^. 



SPIN-2 BEC 



For the case of a spin-2 BEC, the total spin in the two-body interaction F takes on the 
values of 0, 2, or 4, and Eq. ([T]) reduces to [15j 



k k,p,q 



.(2).t 



+^1 ^fmm' ' W'^p,m^q,/i^p+k,m'^q-k,/i' + ^2 ^ ( l)^^^ ttp ,^aq _^ap+k,m^Sq-k,-mMC97) 



where = (4^2 + 3^4)/7, cf ^ = (^4 - ^2)/7, = (7^o - 10^2 + 3^4)/35, and the spin-2 
matrices are given by 



(" 


1 








0^ 
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The order parameter in a spin-2 BEC can be expanded in terms of the spherical harmonics 
of rank 2 as 



1 / 1 ^ 



(99) 



where n =* {nx,ny,nz) and 



Q 



C2 + C-2 - y 3C0 i{C2 - C-2) -m T 1,-1 
t{C2 - C-2) -C2 - C-2 - a/ICo -^(Ci + C-i) 
^ -Ci + C-i "-^^ ' ^ ^ 



-^(Ci + C-i) 

21 



(100) 



2./fCo y 



is a traceless symmetric tensor [151 EH]. In the absence of an external magnetic field, the 
mean-field energy of a spin-2 BEC can be completely characterized by Q, and the ground- 
state phases are ferromagnetic, nematic, or cyclic, depending on the values of the coupling 
constants. The parameter regimes of these phases and the spinor order parameters are given 




2jt 







FIG. 7: (Color online) Phase diagram of spin-2 BECs in the absence of an external magnetic 
field, where C, F, and N show the cyclic, ferromagnetic, and nematic phases, respectively. The 
thick lines represent the phase boundaries. The shape of the wave function is determined from 
\^{9,(j))\^ = I CmY2^(6^, 0)P and the color indicates the phase of ^^(0,0). The axes of the spin 
coordinate are depicted at the bottom right of the phase diagram. Note that the uniaxial (right) 
and biaxial (left) nematic phases are degenerate at the mean-field level. 



by (see Fig. [7|) 

ferromagnetic phase: cf ^ < 0, cf^ > 4cf \ and = (1, 0, 0, 0, 0), (101) 
nematic phase: cf^ < 0, cf^ < 4cf\ and C = ( 0, cosry, 0, ) , (102) 



72 ' ' ' V2 

cyclic phase: cf ^ > 0, cf ^ > 0, and = ^(1, 0, ^2, 0, -1), (103) 

^2") (2") (2") 

where Cq = (45^2 + 35^4) /7, c^ = {qa — 92) H ^ ^2 = (^S'o — 105'2 + 35'4)/35. The nematic phase 
has no magnetization but features a finite spin-singlet pair amplitude, whereas the cyclic 
phase in the absence of an external magnetic field has neither magnetization nor spin-singlet 
amplitude. It is predicted [2^1 [30l [3l] that the ground states of ^^Na is nematic and that 
of ^''Rb lies in the vicinity of the phase boundary between the nematic and cyclic phases. 
(^^Rb is unstable at zero magnetic field in the thermodynamic limit because cf ^ < 0.) The 
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nematic phase has an additional continuous parameter rj (see in Eq. (102)), which is not 
related to the symmetry of the Hamiltonian but represents the degeneracy of the uniaxial 
and biaxial spin nematic phases [SOJETIE^. However, this additional degeneracy arises only 
when the external magnetic field vanishes. In fact, for q < 0^ the biaxial nematic phase has 
a lower ground-state energy than the uniaxial nematic phase, and the boundaries and spinor 
order parameters are given as follows [28j : 

ferromagnetic phase: cf^ > 4cf \ cf ^ < ^, and = (1, 0, 0, 0, 0), (104) 

2n 

biaxial nematic phase: C2^^ < 4cf\ cf'' < — and = ^—^,0,0,0, -> (105) 

(2) . (2) (2) kl (2) 2|g| ^ ,^ /sin0 ^ ^ ^ sin0 

cyclic phase: c\ < Ac) , c\ > — , c\ > , and C = — 0, cost/, 0, ^ 

^ ^ ^ 1 ' 1 2n ^ n ' ^ V 72 V2 

(106) 

mixed phase: > 4cf \ cf ^ > ^, > and = (cos^^ 0, 0, sine^ 0) 

2n n 

or (0,sine^O,0,coseO. (107) 



where cos0 = y 1/2 + qjnS^^ and cos^' = — q/?)nc{ \ The phase diagram is shown 

in Fig. |8](a). The cyclic phase for g < does not have any magnetization but has a finite 
spin-singlet pair amplitude that depends on q. The mixed phase does not have a spin-singlet 
pair amplitude but has a finite longitudinal magnetization that depends on q. 

On the other hand, for g > 0, based on a numerical calculation, we find that the ground- 
state phases are given as follows: 

broken-axisymmetry phase: = (±a, 6, c, 6, ±a), + (— ) sign for c^^ < (> 0) (108) 

uniaxial nematic phase: = (0, 0, 1, 0, 0), (109) 

11 .f^f^iiiO ^ sin^\ 

cyclic phase: Q = I — ^,0,cosc/,0, ^ J , (110) 

V v2 v2 / 

(2) (2) 

where a, 6, and c are positive except for the case of = and C2 > 0, where a = 0. The 
broken-axisymmetry phase has a transverse magnetization and a spin-singlet pair amplitude, 
both of which depend on g, as in the case of the spin-1 broken-axisymmetry phase. The 
phase boundaries between the broken-axisymmetry and cyclic phases and those between the 
cyclic and uniaxial nematic phases correspond to those between the mixed and cyclic phases 
and between the cyclic and biaxial nematic phases, respectively. The boundary between 
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the broken-axisymmetry and uniaxial nematic phases is determined numericaUy. The phase 
diagram is shown in Fig. [sF 



{a)q<0 

C2 n 




FIG. 8: (Color online) Phase diagrams of spin-2 BECs for nonzero q. (a) Case of g < 0, where BN, 
C, F, and M stand for the biaxial nematic, cyclic, ferromagnetic, and mixed phases, respectively. 
In the absence of an external magnetic field = 0), the C and M states can be transformed to 
each other by a rotation in space, (b) Case of g > 0, where the UN and BA stand for the uniaxial 
nematic and broken-axisymmetry phases, respectively. The thick lines show the phase boundaries. 



The Bogoliubov Hamiltonian in a spin-2 BEC is given by [15j 
H.^ = ^(cf^ + cf^f)^ + 4cf)|(._)p)+,iV((rn 

(ek - ncf\if - 4ncf K._)p + qm' - q{{r?)) 



nc, 



(2) 



(2) 

° ' 2DId^ + I)kl)-k + ) + ^ ( 2Fi, • Fk + Fk • F_k + fL • F^.^ 



2 V - ■ K -Kj ■ 2 

+ncf\f) • fmm'alc,,„ak,m' + 2nc^2 \S^)^S^ 

+ncf\-ir ((s-)«U«-k,-m + {s-r&Kma-K-m) 



(111) 



where 



m 

is the spin-singlet pair amplitude, and 



(112) 



(113) 
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is the pair fluctuation operator of the condensate. As in the case of a spin-1 BEC, we 
substitute for c>^^ {% = 0, 1, 2) in the sum over the momentum in Eq. (111). 



A. Ferromagnetic phase 



For the ferromagnetic phase (101), the sub-Hamiltonian with respect to ak,2 is diagonal- 



ized by the Bogohubov transformation [L5\ 



Ok,2 = \l 7^ «k,2 + ■ 



.(2) ^ , J2), 



2-E'k,2 

with the Bogohubov spectrum 



2E, 



k,2 

—a' 



k,2 



k,2' 



k,2 



fk 



ek + 2n(4'^+4cf^) 



(114) 



(115) 



The diagonahzed Hamiltonian is given by 



^eff 



-E'k,2fek,2^k,2 + (ck - 3g)a^ak,i + (ck - 4g - 4ncf ^) ak,o^k,o 



+ (^fk - - 6ncf^ j o-k -i^k -1 + (ek - 8ncf ^ + 2nc2^^^ aj^ _2«k -2 



(116) 



where 



< = 4qN + 



:(2)> 



5]'(6k + n(4')+4cf))-^^k,2) (117) 



is the GSE. The m = 2 mode is massless, and in the absence of an external magnetic fleld, 
the m = 1 mode also becomes massless. For the excitation energies to be positive, we must 
have + 4cf^ > 0, g < 0, cf^ < |g|/(2n), and > 4cf\ The GSE , pressure P, 
sound velocity c, and quantum depletion (A^ — No)/N can be calculated as 



El 
V 



M 



1 + 



128 na 



15 V TT 



27rh^n^a4 
M 



M2 



1 + 



1 + 



64 /na^ 



TT 



TT 



N - No 8 /na 



TV 3 V TT 

where the renormalization of is carried out to obtain the finite GSE. 



(118) 
(119) 
(120) 
(121) 
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B. Nematic phase 



Let us next discuss the nematic phase. We discuss the cases of q = and q separately, 
since the stationary solutions are different. 



1. Case of q = 



In the absence of an external magnetic field, the nematic phase is characterized by Eq. 



(102). We introduce the following four independent fiuctuation operators: 

sinry 



" V2 



(«k,2 + Ok -2) + cos?7ak,o, 



V2 



^Kfy — yl*^"^^'^ ~^ -1), 

where d\^^d and CL\^j^^y^^ describe density and spin fiuctuations, respectively 
we introduce the following nematic fiuctuation operator: 

cos 77 



-(ak,2 + ^k-2) -sinr7ak,o. 



(122) 

(123) 

(124) 

(125) 
Furthermore, 

(126) 



This operator describes fiuctuations of the nematic order parameter; in fact, the coefiicients 



(127) 



in Eq. (126) are related to the nematic order parameter (102) by 

f COST] COST] 

In terms of the operators in Eqs. ( 122[ )-( [T26 ), the total Hamiltonian can be decomposed 
into sub-Hamiltonians as 



H ^ l«k,d^-k,d + «k,d«-k,dj 



+ 
+ 

+ 



(2) 



ei, + ncf\r] - n/3) 



(?7 + 7r/3)] aljjkj. 



ncl^V?7 + 7r/3) -^t . . > 
2 (</.«-k,/. +«k,/.a_k,/., 

ncl^V?7 - 7r/3) , i * . . , 
2 (</,«-kj, + «k,/,a_k,/J 



+ 



(2') 

ek + n4 (??) 



.(2) 



k,/,"'-k,/^ 



ak,/,a-k,/J 



(2) "I 

/ (2)\ . ncn - - \ \ 

+ (ek - n4 'j <^ak,„ + ^-«,^«-k,r, + «k,„a-k,^) >, 



(128) 
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where c"^\r)) = 40^"* sin^ r) — To diagonalize this Hamiltonian, we consider the foUowing 
Bogohubov transformations: 



(2) , J2)^ 



£^1 



-«k,d + 



a' 



k,d 



k,(i' 



ek + nci\r) + 7r/3) + E^^j^ ^ , , / ^k + ^4^^ (v + 7r/3) - £;k,f. 



(2)/ 



5k,/. 



£^1 



-«k,/. + 



£^1 



k,/x 



ek + ncf (77-7r/3) + £ ;k,/, 



-«k,/. 



-k,/y' 



ek + n4^^(?7) + £;k,. 



ek + ncf^(?7) -£;k,/..t 



£^1 



-«k,/. + 



E, 



ek - nc'^^ + ^^k,,; . / ek - nc^^^ - ^^k,,; 



(2) 



E, 



k,?7 



-«k,)? 



£;i 



a' 



k,?7 



k,J7' 



where the Bogohubov spectra are given by[20l [2T] 

Eu,d = \/ek[ek + 2n(4')+cf )], 



£^k,/. = Vek[ek + 2ncf (77 + 7r/3)], 
^k,/, = V'ek[ek + 2ncf (77-7r/3)], 
^k,/. = V'ek[ek + 2ncf (?7)], 
£^k,,; = V^k(ek - 2ncf^). 



(129) 
(130) 
(131) 
(132) 
(133) 

(134) 
(135) 
(136) 
(137) 
(138) 



The diagonalized Hamiltonian is then given by 



Eu,dblJ>k,d + £^k,/Ji,/Jk,/. + £^k,/,^i,/^^k,/, + £^k,/,^i,/,^k,/, + £^k,„^lc,^^k,„ 

(139) 



where 



E 



(ek + n{cf + ) - £;k,d) + (ck + ncf{r] + 7r/3) - Ek,. 



+ (ek + ncfij] - T^/?>) - Ei,j^^ + (ek + ncf (ry) - £;k,/.) + (ek - ^4'^ - £^k,r,) 



(140) 



Although the total number of the symmetry generators of the Hamiltonian is 4, the five 
Bogohubov massless modes emerge. The physical origin of this result is discussed in Sec. 
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V. We also note that regardless of the value of 77, the above Bogoliubov spectra are real 
if Cq^^ + c^2^ > 0, C2^^ < 0, and c^"^ < 4c^^\ in agreement with the stability criteria of the 
mean-field ground state. As shown in Appendix A, the GSE E^{r])^ pressure, sound velocity, 
quantum depletion are respectively given by 



V 



«'(cr + cr) ( 16VM3 



+ [n(2cf) - )]i + Xcos(277 + 27rj/3)]l \, 
i=o 



(141) 



1 + 



8VAP 



f)3 



+ 



4VM3 



(n|cf I 



+ [n(2cf) - )]i ^[1 + Xcos(2ry + 2nj/3)f^ \, 



j=0 



(142) 



C{7]) = 



'n(c?^+cg 



1 + 



v/n(cf ) + )3 



n(- 



,(2) , ^(2) 



-0 



+ ' J ^(1 + X cos(2r7 + 27rj/3))i 



,(2) , J2) 



-0 



+ 4 



j=0 



N-No 
TV 



VM3 
3^ 



Vn(c?) + )3 + (cf (r; + 7r/3))' + {cf\v - tt/S))' 

+ ^n(cf) (r7))' + Vn(|cf|)3 



(143) 



, (144) 



where X = — 2cf^/(2cp^ — 4^^) and 2cf^ — > in the nematic phase. 
We note that the GSE satisfies the following symmetries: 

E^iv + ^/3) = ^o^(^) and E^{-r)) = E^{r)). 



(145) 



Therefore, other physical properties derived from it such as the pressure, sound velocity, 
and quantum depletion also satisfy similar relations: we may therefore restrict the domain 
of the definition of to < r; < 7r/3 to find the ground state. We define the part of the 
GSE that depends on rj as 



E/cj = J^[l + Xcos(2r^ + 27rj73)]t, 

j=0 



(146) 
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where uj = ^\/M^V[n{2c{^ /(157r^/z^). Since = and 7r/6 satisfy dE/dr] = 0, these 

points are the candidates of the ground state. The second derivative with respect to rj is 

1 d^E 
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FIG. 9: Second derivative of E with respect to 77. The solid and dashed curves correspond to 
77 = and 7r/6, respectively. 



plotted in FIG. [9| which shows that = and = 7r/6 correspond to the ground state for 
X < and X > 0, respectively pOl [2T] . The phase transition between uniaxial and biaxial 



nematic phases occurs at X = because the diagonalized order-parameter matrices (100) 
for X < and X > given respectively by 



Q 



X<0 



V 





















O /2 . 



and Qx>o 





^2 



(147) 



cannot be connected by any transformation that belongs to x aS0(3) because 

Det((5x>o) = and Det((5x<o) 7^ 0. The phase diagram that incorporates the quantum 
fluctuations arising from the Bogoliubov theory is depicted in Fig. [TOj Here, we note that 
the many-body analysis in the single-mode approximation (SMA) [151 ESI ES] does not pre- 
dict such a phase transition because the SMA ignores quantum fluctuations of momentum 
that play a crucial role in the phase transition between the uniaxial and biaxial nematic 
phases. 

We briefly comment on another stationary state = ^^(0, 1, 0, 1, 0). This state does not 

v2 

provide any new piece of information because it can be transformed into a biaxial nematic 
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FIG. 10: (Color online) Beyond-mean-field-theory phase diagram of spin-2 BECs in the absence 
of an external magnetic field obtained by the Bogoliubov theory, where C, F, UN, and BN stand 
for the cyclic, ferromagnetic, uniaxial nematic, and biaxial nematic phases, respectively. The thick 
lines show the phase boundaries. The phase diagram is consistent with that obtained in Refs. 

[201 EI]. 



state by an element of U{1) x S0{3)] in fact, 

C = e^^e-^^^e-^^^'e-^^^^C^^, (148) 

where = ^(1, 0, 0, 0, 1), 9 = 7r/2, a = 7r/4, /3 = 7r/2, and 7 = 0. Therefore, this state 
v2 

belongs to the biaxial nematic phase and it has the same spectra as the biaxial nematic 
state in the absence of an external magnetic field. 

There is a hidden symmetry in the mean-field stationary solution of the nematic phase. 
The nematic phase is characterized as \{s-)\ = 1/2 and (f) = 0, both of which remain 
invariant under the U{1) x S0{3) transformations and are independent on rj. In fact, we 
can show that the full symmetry group that leaves the nematic invariant is U{1) x 5^0(5), 
which includes the U{1) x 5^0(3) as a subgroup. This is confirmed by using the fact that 

_(2) 

the pair-singlet interaction term, which is proportional to C2 in the Hamiltonian, has the 
5^0(5) symmetry [33j, and that (f) = if \{s-) \ = 1/2 [ISll^J. In Sec. VI, we discuss a role 
of the SO (5) symmetry and relationship with the massless modes in the nematic phase. 
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2. Biaxial nematic phase 



The effective Hamiltonian of the biaxial nematic phase defined in Eq. (105) can be 



decomposed by the transformations of Eqs. (122)- (126): 



ttBN 
^eff = 



2 

.(2) ^ J2)x 



H [Cik4^-k4 + %,d<^-k,^i; 



.(2) J2)x 



j — fx,fy 



He, + 4nc?-nc?)al,a,j^ + 



2 ■("kj.^-kj. +4,/,«-k,/J 
+(ek - 4g - n4^0«k,„«k,r, + ^^(ak,^a_k,^ + 4,„«-k,^)|- (149) 
This Hamiltonian can be diagonalized by the following Bogoliubov transformations 

bk,d = 



-.(2)\At 



nc 



(2) 



-«k,d + ■ 



£^1 



^1 



k,d . t 

—a' 



k,d 



k,(i' 



(150) 



. . / ek - 3q + n(cf^ - cf ) + E^j^ , / ek - 3g + n(cS^^ - ) - ^k,/, n^i^ 
Ok,i = \l ^— akj + y ^r^^ a_^^., (i5ij 



kjt 



£^1 



kjt 



'ek + n(4cf)-cf) + Ekj. 



6k + n(4cf)-cf)-gk,/..t 



£^1 



k,/z 



■«k,/. + 



£^1 



kj. 



-k,/.' 



^k,f; 

with the result 



' fk - 4(? - n4^^ + ^^k,,; 



'ek-4g-n4^^ -^k.r/.f 



£^1 



Q-k,! 



kfl 



£^1 



k,»; 



-k,7?' 



(152) 



(153) 



= ^o"""^ + E [^k,ii,/k,d + E^jXhlfMh + 4,/,^k,/J + EkjAjXf. + ^k,,&i,,&k,,] , 

k 

(154) 



where 



E 



BN 



^ek + nco^^ + nc2^ - E'k,, 



-2 (ck - 3g + ncf ) - ncf - Ek,/,) 

- (ck + 4ncf ^ - ncf - Ek,/,) + (e^ - 4q - ncf^ - Ek,^ 



(155) 
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is the GSE, and the Bogohubov spectra are given by 

Eu,d = 



ek(ek + 2n(cf +cf )), 
£^k,/. = V (^k - 3<?)(ek - 3g + 2n(cS') - )) 



ek(ek + 2n{4c\ 



(2) ^(2) 



)), 



;ek-4(/)(ek-4(/-2ncf ). 



(156) 
(157) 
(158) 
(159) 



The density and spin fluctuations around the z axis are massless and those of nematic and 
spin fluctuations around the x and y axes are massive. This is because the symmetry of 
the Hamiltonian is reduced to U{1) x S0{2) due to the external magnetic field and the fact 
that the isotropy group of the biaxial nematic phase does not include any continuous group. 



The Bogoliubov spectra in Eqs. (157) and (159) are positive semideflnite if g < 0, which is 



consistent with the stability criterion of the mean-field ground state. 

The GSE Eq^ , pressure, sound velocity, and quantum depletion up to the LHY correc- 
tions are given as follows: 



ET 
V 



^n{cf^ + )3 + ^^,{{n\c?YHi{h - sgn(cf )) 



+2(n|cf) - 4'^|)i0i(t4 + sgn(cf) - )) + [n(4cf) - )]§}, 



(160) 



n2(4^^ + 4' 



,(2)> 



1 + 



8VM3 



,(2)^ 



+ ^^{(n|cf |)t [301 (ts - sgn(crO) " 2(^5 + iWAh - sgn{c'^>)) 
+2in\cf^ - |)t [SMU + sgn(cf) - )) - 2{U + + sgn(cf) - c^^) 

+3[n(4cf)-cf)]t}, 



(161) 



'^^^°'+^'^)Jl + ^;/n(4ncr))3 + 



M 



/M3 \Jn\cf\^ 

157r2^3 (2) (2) 
cq -r C2 



'/>2(t5,-sgn(cf^)) 



+ 



2VM3 V^Ih 



,(2) J2)| 



157r2;i3 c^2) (2) 



(.2(t4,sgn(cf)-cf)) 



TT 



cq -r ^2 



(162) 



N 



VM3 
3^ 



^(e(2) + )3 + 2Vn|4^^ - c^f + n{\^^'\fUh - sgn(crO) 



.(2) _ J2)|3 



,(2), 



,(2)^ 



+03(^4 + sgn(ci') - cf)) + Vn(4cr - 



.(2) _ J2)x3 



(163) 
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where = —3q/{n\cf'^ — ^^^1) — 1 and = — 4g/(n|c2^^|) — 1. For g < 0, all the LHY cor- 
rections are positive definite, indicating the robustness of the biaxial nematic phase against 
the beyond-Bogoliubov eflFect. 

3. Uniaxial nematic phase 



The effective Hamiltonian of the uniaxial nematic state in Eq. (109) is obtained by using 



the canonical transformations of Eqs. (122)- (126) as 



12. 



UN 
eff 



\ (ek + ncy + nc\ 0«M«k,d + ^ (ak,da-k,d + <d«-k,d) 

n(3cf)-cf; 



|^(ek + g + 3ncf ^ - nSf^)a}^ -a\^^j + 



j — fx,fy 



+ + 4:q- ncf )a[ ;ak,i + -^{di^,id-u,i + «lc,/«-k,i) 



(2) 



(164) 



This can be diagonalized by means of the foUowing Bogohubov transformations, 



, /ek + n(4')+cf ) + ^k,d. , /ek + n(4'^+4'0-^M.t 
Ok,d = \ «k,d + \l « 



.(2) , J2)x 



k,d 



k,d' 



(165) 



, ; ek + g + n(3cf)-cf ) + Ek,/ / ek + g + n(3cf ^ - ) - ^k,/. nfifi^ 



^1 



&k,Z 

with the result 



E, 



Cki 



k,?7 



£^1 



k,77 



(167) 



ht = Er + E [^M%A^ + ^k,/.(&i,;jk,/. + %/,^kjj + ^M(&i,;ikj. + &i,,&k,,)] , 

k 

(168) 



where 



E, 



UN 



(ek + n4^^ + ncf ^ - ^^k.d) 



-2 (ck + <? + 3n4') - ncf - Ek,/,) + 2 (ck + 4g - ncf - £;k,^) 



(169) 
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is the GSE, and the Bogohubov spectra are given by 

Eu,d = \/ek(ek + 2n(4') + cf )), 



£^k,/. = V (^k + <?)(ek + q + 2n(3cS') - )), 
Ek,^ = \/(ek + 4g)(ek + 4<7-2ncf ). 



(170) 
(171) 
(172) 



The density fluctuation is massless, while the nematic and spin fluctuations around aU the 
axes are massive. This is because the isotropy group of the uniaxial nematic phase includes 
the continuous group S0{2) and, therefore, one massless mode should appear. As can be 
seen from Eqs. ( 17lD and ( |172 ), the Bogoliubov spectra are positive deflnite only if g > 0, 
which is consistent with the mean-fleld analysis. 

The GSE Eq^ ^ pressure, sound velocity, and quantum depletion are given by 
UN ^2^^(2) , ^(2)^ / ^a.Hwrs i — \ 



+{n\3c? - |)Mi(^6 + sgn(3cr - c^)) 



(2) ^2)^ 



(n|4^^|)t<^i(t7-sgn(4^^)) 
(173) 



n2(4^)+crO 



P = 



1 + 



8VAP 
5^ 



(2) 




+ 0^)3 



+^^{(n|cf 1)1 [301 (tr - sgn(4^))) - 2(t, + lW,{t, - sgn(cf ))' 
+(n|3cf) - cf |)i [301 (te + sgn(3cf) - )) - 2{U + l)0'i(t6 + sgn(3cf) - ))] }, 

(174) 



c = 



n (4ncf) 
M 



n\c 



(2) I 



157r2^3 \ ^(2) ^ J2) 



02(^7, -sgn(4^^)) 



n|3q 



(2) J2)| 



J2) (2) 
-0 "T ^2 



-02(^6, sgn(3cf)-cf)) 



(175) 



7V-7Vo 



VM3 
3^ 



^(c(2) + c(2))3 + 2Vn(|cf |)303(t7 - sgn(cf )) 
+2X(|3cf)-c?)|)303(t6 + sgn(3cf ) - )) 



(176) 



where te = 1 3cf ^ — C2^^ |) — 1, ^7 = 0^/(^1 ^2^^ | ) — 1 • The LH Y corrections become imaginary 
for g < 0, which implies that the system undergoes the dynamical instability, and that the 
uniaxial nematic phase is only stable for g > 0. 
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C. Cyclic phase 



In addition to the stable configuration of Eq. (106), we also examine the tetrahedral 



configuration of Eq. (103), even though the latter configuration is not a stationary solution 
of the mean-field theory for nonzero q. This is because the tetrahedral configuration has 
attracted considerable attention, since it gives rise to nontrivial phenomena such as the 
non-Abelian vortices 



1. Stable configuration for nonzero q 



For the stable cyclic phase in Eq. (106), we consider the following canonical transforma- 
tions: 



sin^ 



(ak,2 - + cos6>ak,o, 



V2 

sm COS0 

sin + 6 cos^ 6 

— sm cos^ . 

I , ^ — Ok I — i 

V2 sin^6' + 6cos2^ ' 

1 . , 

^K,2 + ak,-2), 



— sm COS 9 , 



V2 sin^ e + 6 cos^ 9 
sm ^ + ^3 cos^ 
V2 sin^ ^ + 6 cos^ 9 



cos 9 

71" 



(-«k,2 + Ok -2) - cos6'ak,o, 



(177) 

(178) 

(179) 
(180) 
(181) 
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where dk^a, (^kj^, ^Kfy^ ^^'^ ^Kfz represent the density and spin fluctuations, and dk^e is the 
operator that commute with the other four operators. Then, the Hamiltonian is given by 



= 2qN + 



ek - g + 2ncf ^ ( 1 + 



nc. 



(2) 



+ 



ncf^ I 1 



nc. 



(2) 



sin^^ + Scos^^ 



sin2^ + 3cos2^ 



f Uj 

jy ^-iJy^ 



2V3gsin26' . ^\ ^ I ^ (2) ^ ^qcf^ \ .t 

(ak,/.a-k,/, - </,a-k,/J + Uk + 2nci > ^2q ^ </,«k,/. 



sin^^ + 3cos2^ 



1 



+- I 2ncf^ + 2q 



(2) 



T(2)" 



^(^1. -I- 75./^?^ — n.A-^^ 



(ak,/,a_k,/, +at_^^aLk,/J 



(ok,da-k,d + ^k,d^k,d) 



-2qsm20{dl.^dk^0 + dlgdk^e + dk,dd^k,e + ak,d^Lk,e) }, 



(182) 



where Cq^ = 4q'^/{n'^c^^). For nonzero q, the above Hamiltonian is not decomposed into 
sub-Hamiltonians completely because all modes except for d^j^ couple. For the mode, 
we consider the following Bogoliubov transformation: 



ek + 2ncf^ + 2q- 4:c\'>q/4'> + E^j. 



(2)^/J2) 



^kj, = 



ek + 2ncf ) + 2(7 - 44^^g/4^^ - E^j^ 



(2)^/J2) 



2E^ 



2E, 



where 



y^k (ek + 4ncf ) +4q- Scf^g/cf ^ 



(184) 



is the Bogoliubov spectrum. For the fx and /^^ modes, we consider the following Bogoliubov 
transformations: 



^k,/t] 



1 



'ek-g + 2ncf) + 2cf)g/4') 



£^1 



£^1 



(ak,/, +ak,/J 



'ek - g + 2nc? + 2cf ^g/cf - ^k,/< 



-(^^-kj.+^lk,/,) 



£^1 



k,/t 



(185) 
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(186) 



nc. 



sin ^ + 3 cos^ sin 9 cos ^ 

sin^ ^ + 3 cos2 e 



(187) 



where r = — Im/3/Re/? and 

£;k,/, = ^l^n^cf^cf^ - 3q^ + ek (ek + ^ncf^ + Ancf\/c^p - 2q) , (188) 

is the doubly degenerate Bogohubov spectrum. FinaUy, for the density and 9 modes, we 
consider the foUowing transformations: 



where Bk =* (5k-,&k,+), ^^do - («k,d, ak,(?) 
U{k) 



if 2cU)E.ik) +C,-{k)X,^{k) 4gekC.-(fe)sin2^- 3^^^^Jg^;j,),,„,, 



2C2+(fc)g2(fc) + <^2+(fe)X2-(fc) 4?ekC2+(fc)sin26' + 



(189) 



(190) 



V{k)^ 



+ C2- {k)X2+{k) 4(7ekC2- ( fc) sin 2^ + 



X2-(k) 



\ I 2C2_(A:)E2(A:) ' ^^-V'V^»-i:+V'V ^'i>-K^z- w '^■^^^ ' 8gekC2- (fc)S2(fc) sin 26» 

2^ 2cOT»)+^2+(fc)^2-(/c) 4g6kC2+(/c)sin2^- 3^^^^Jffg;,)^,„,, 



^2±(A;) = -n(4') - 2cf + 2cf ))ek + 2n2cf (cf - ) ± E^ik), 



C2±{k) 



E 



'k,± 



^2%(^)ek + 16g2 sin^ 2^e^(ek + 2nci^) ' 



(191) 

(192) 
(193) 



E^ik) = {(ncf ) - 2ncr )^ + 4n^cf (cf - c?))}^^ - An'c^\c^^ - 2c?))(c?) - cf ))ek 



1/2 



+{2n^cr(cr-cf))}n (;L94) 



and £^k,d= are the Bogohubov spectra given by 



^k,± = ^ek (ek + n4'^ + 2ncf ) + 2n^cf^ (4^^ - 4'^) ± E, 



(k). 



(195) 
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Using the above transformations, the total Hamiltonian is diagonalized as follows: 



where 



k 



2 (ck + 2ncf) + 2cf -q- E^^j^j 



+ (ek + 2ncf ) + 2g - Ac^^q/cf - E^j}j + (2€k + ncf + 2ncf - ^k,+ - ^k,-) 

(197) 

In the long-wavelength hmit, E'^ _ c^i 2nc^^^e\^] therefore, - is hnear and massless. This is 
the NG mode associated with the spontaneous breaking of the U{1) gauge symmetry. The 
spectrum on the spin fluctuations around the z axis is also linear and massless since the 
SO {2) symmetry is spontaneously broken at the mean- field level. The spectra around the 
transverse axes are massive because the rotational symmetry about transverse directions are 
explicitly broken by an external magnetic field. We note that this cyclic configuration is 
stable if Cq^^ > < 4cf\ cf^ > \q\/{2n)^ and > 2\q\/n. It is robust regardless 

the sign of which is consistent with the stability condition of the mean-field ground state. 
Using the renormalization procedure discussed in Appendix A, we find that the GSE per 
volume V is given by 



rf = 2^^+^r k + 



^2 



(2ncf ) + 2q- 4cf)g/cf )i + 2(2ncf ))i07±(t8) 



+ (n4'))t 08(^9) 

where the + (— ) sign in 4'7±{t) corresponds to the case of g > (g < 0), 

4±it)-Z2{t) 



,(198) 



15 { 

0v±(t) = y dxx^ + 2i±(t) - v^^2(^)+a;2(a;2 + 2^i±(t)) 



2x^ 



(199) 
(200) 



0?^(t) = )Jx^ + Z3X^ + z,{t) ± ^Jze{t)x^ + Z7{t)x^ + zl^), (201) 
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with ts = \q\/2ncf\ t, = {2qf/{ncf)\ z,±{t) = 1 Tt(l - 2cf V^f ), ^2(t) = t\l2cf^/c^^^- 
3), zs = (4^^ + 2cf z,{t) = t{l + 4^Vcf ), z,it) = 2(cf - 2t, z,{t) = (1 - 

2cf /c(^))2 + 4t(4^Vcf - 1), and z^t) = -4(cf /c(^))(4^) - 2cf )[cf /(ci^))^ - t/cf ]. The 



plots of (/>7±(t) and 08 (i) in Fig. |ll(a) show that the spin components of the GSE around 
the X and y axes for g < (for q > 0) decrease (increase), while the density component of 
the GSE increases as the external magnetic field increases. The typical values of ts and tg 
are of the order of 10~^ and 10~^, respectively, for the parameters of ^''Rb. Performing the 




0.01 



0.02 



0.03 



0j04 



l.xlO 2.x 10 3.x 10 ' 4.x 10 ' 



FIG. 11: Plots of (a) (f)7±{t) and (b) (f)s{t)- These functions are approximated as (f)7±{t) 2^ 1 zb 
24.lt + 116^2 and Mt) - 1 + 1140t. 




FIG. 12: Plots of (/)9±(t) (a) and (/)io(t) (b). 



derivatives of the GSE with respect to V and n, we obtain the pressure and sound velocity 
as follows: 



2 (2) 



2q^ 4VM3 
2 ^^^15^ 



{2ncf'' + 2q - 4cY>q/c\">)^ 



nci + q — 2ci q/c2 



(2)^/J2), 



(202) 
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(a) 




0.01 0.02 0,03 0*4 



1.x 10-* 2.X10-' 3.x 10"* 4.x 10-' 



FIG. 13: Plots of (a) (f)n±{t) and (b) (f)i2{t)- These functions are approximated as (/)ii_(t) 
1 - 15.lt - 40^2, (/>ii+(t) ~ 1 + 13.6t - 26*2, and (puit) ~ 1 + 700t. 



' nc, 



(2) 



where 



M 



1 + 



4VM3 c\ 



(2) 



n{2ncf^)^,Ats) 



4VM3 c^- 



(2) 



nc 



(2) 



7^2/^3 4'^ 2ncf^ +2q- Ac^^^q/c 



(2) 
2 



(203) 



09±(t) = 1507±(t) - 12t0'7±(t) + 4t^0'7'±(t) 



(204) 



(/>io(t) = 1508(t) - 16t(/>'8(t) + mmit). 



(205) 



The behaviors of (/>9±(t) and (/>io(t) plotted in Fig. 12 show that for g < the sound velocities 
with respect to the density and spin wave around the x and y axes decrease, while for g > 
those with respect to the spin around the x and y axes increases, as the quadratic Zeeman 
effect becomes stronger. The quantum depletion is expressed as 



n ( 2q/n + 2cf ^ - 4cf )g/(ncf ) ) + 2 Vn(24^0Vii±(t8) 



.(2)^3, 



+ V'n(4'))30i2(t9) 



, (206) 



where the + (— ) sign in corresponds to the case of g > (g < 0), and 



(l)ii±{t) = / dxx'^ 



X 



lTt(l-2cfVcf) 



3t2(4cfVcf - 1) + x^{x^ + 2 T 2t(l - 24'74'0) 



(2)/J2)^ 



^12 



(t) = ^ rf^^'(^n + ^il + ^2^1 + ^2l)- 



1 ,(207) 



(208) 
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The behaviors of (/>ii±(t) and (/>i2(t) plotted in Fig. 13 show that the quantum depletion from 
the density fluctuations increases, while that from the spin fluctuations around the x and 
y axes for g < (for g > 0) decreases (increases), as the external magnetic fleld increases. 
The variations of the quantum corrections with respect to the density component are small 
since the changes in 08 (t), 0io(^)5 and 012 (i) are of the order of 10~^. On the other hand, 
the variations with respect to the spin components are of the order of 1. Considering the 

(2) ('2) ('2) 

fact that Cq ^ c\ \ C2 , and \q\/n for the alkali species, however, the main contribution in 
the LHY corrections stems from the density fluctuations as in the cases of the other phases. 

2. Tetrahedral configuration 




FIG. 14: (Color online) Schematic illustration showing that the symmetry axes of the cyclic order 
parameter constitute a tetrahedron, where the color shows the phase of the order parameter as 
indicated in the color gauge of Fig. 7. 



Finally, we examine the mean-fleld conflguration of Eq. (103), which has the tetrahedral 
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symmetry as shown in Fig. [T4| We consider the foUowing canonical transformations: 



«M = ^(«k,2 - «k-2) + ;^ak,o, (209) 

" 2^ «k,-i, (210) 

akj, = ^ 2^ «k,i-^ ak,-i, (211) 

«k,/. = ^(«k,2 + ak,-2), (212) 

flk.p = ^(-«k,2 + «k -2) + ;^^k,o, (213) 

which represent the density (d), spin {fx,fy,fz), and pair (p) fluctuations, respectively. 
Then, the Hamiltonian is rewritten as follows: 

.2^(2) 



= 2qN + + Yl [ (^k - ^ + 2?^cf ^) 4,i«k,j + ?^cf ^(«k,,«-k,, + <i«tkj) 

k I j=fx,fy 

+ (ek + 2g + 2cf^) 4,/,«k,/. + ^cf^(4,/,«-k,/, + «k,/.«-k,/J 

+ (^Ck + nc[)^^) «k,d«k,d + (ck + 2n4^^) «k,p«k,p - 2q{al^^ai^^p + a^pak.d) 



(2) 

+^(4,>-k,d + «k,da_k,d) }>. (214) 



For nonzero the density fluctuation couples with the pair fluctuation. For the spin modes, 
we consider the foUowing transformations: 



ek + 2g + 2ncf ) + E^j^ + 2q + 2ncf^ - E^j^ . 
= Y «k,/. + y ^^^^^ a_^j^, (216) 

where the Bogohubov spectra are given by 



^k,/. = y (ek - g) (fk - + 4ncf , (217) 
^kj. = ]J{eu + 2q) (eu + 2q + 4ncf^y (218) 
On the other hand, for the density and pair modes, we consider the foUowing transformations: 

Bk = [/(A:)Ak,dp + l^(fc)At_k,ap, (219) 
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where Bk =* (6k-,fek,+), Ak,dp =* (ak,d,«k,p), 

+ C3-(fc)X3+(A:) 4C3-(fc)(7(ek + ncf ) - 



X3-(fc) 



C(fc) = 2 



2C3-(fc)E3(fc) 



- 2C3+(feV3W +^3+(fc)X3-(fc) 4C3+(fc)g(ek + ncf) 



■7jrr-rurw-nT + C3-(k)X3:(k) 4Cs-(k)q(eu + ncf^) H 



SC3+(k)q(eu+nc^2^)E3{k) 



Xs-jk) 



Vik) = - 



2CsJ)Es(k) + Cs+ik)Xs-{k) 4C73+(fc)g(ek + nc','^) 



8C3+(k)q{e^+ncf'>)E3(k) 



X3±{k) = -n(4') - 2cf )ek + 2(ncf )2 ± E^ik), (222) 
Esik) = y {n(c(') - 2cf )ek - 2(ncf )2}' + Wq^e^ + nc^^^) [e^ + n{c^^^ + )}, (223) 

C3±(fc) 



.(2) 



(2)n2 



£^1 



k,± 



\\ Xi^{k)e^ + 16g2(ek + ncf ) |(ek + 2ncf )(ek + ncf ) - X3^(A;)} 
where the Bogohubov spectra are given by 



(224) 



E^,± = \lel + n{cf + 2cf )ek + 2{nc'^'^f + Aq^ ± £;3(A;). 



(2)n2 



(225) 



By using the above transformations, the effective Hamiltonian is diagonahzed as follows: 

^eff = + e' [^kj.(&ijjk,/. + bij.Kfy) + ^kj.&Lj, V/. + ^k,+&i,+6k,+ + ^k, Ji, Jk,-] , 

k 

(226) 



where 

EC = 2qN+—^--Y, 



2 (ek - (? + 2ncS') - E^,j)j + (ek + 2(7 + 2ncf ^ - £;kj.) 



(2ek + ncj,^^ + 2n4^^ - £;k,+ - E^, 



(227) 



is the GSE. As can be seen from Eqs. (217) and (218), the cyclic configuration of Eq. (103) 



suffers the dynamical instability unless q = Moreover, the following inequality needs to 



be satisfied for the Bogohubov mode of Eq. (225) to be stable: 



.(2) A. 3 



2j2) J2) 



.(2)^2 



,21,2 



+ [^it'cficff - 16g^n(4^^ + 24^0]ek + m\q' - n'c)^>c\'>) > 0. (228) 



.(2) 



,(2) 



.2/^2 ^2 (2)J2)^ 
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Conversely, the dynamical instability occurs if the above inequality is not satisfied. There- 



fore, the tetrahedral configuration (103) is unstable for an infinitesimal quadratic Zeeman 



effect in the thermodynamic limit. For q = 0, the spectra (217), (218), and (225) reduce to 
those of Ref. [E]: 



Ekji = y Ck (^Ck + 4ncf , {i = x, y or z) 



ek + 2nc2^^ . 



(229) 

(230) 
(231) 



The above spectra are positive semidefinite if Cq > 0, ^ > 0, and C2 > 0, consistent with 
the stabihty criteria of the mean-field cychc state. 



VI. RELATION BETWEEN THE NUMBER OF THE SPONTANEOUSLY BRO- 
KEN GENERATORS AND THAT OF THE NAMBU-GOLDSTONE MODES 

We discuss the number of the Nambu-Goldstone (NG) modes in hght of the rule found 
by Nielsen and Chadha [22j (see Ref. [34] for a review). In general, if the Hamiltonian 
has a symmetry with respect to the internal degrees of freedom, we can introduce the 
corresponding conserved charge operator defined as follows [35] : 

G = J d^j\x), (232) 

where 

jO(^) = l^^^sif^ix), (233) 

S is the action corresponding to the Hamiltonian, 5^^^(x) is an infinitesimal transformation 
of the field, and the overdot denotes the differentiation with respect to time. The conserved 
charge operator commutes with the Hamiltonian, 

[G, H] = 0. (234) 

When the spontaneous symmetry breaking occurs, the ground state |GS) does not have the 
full symmetry of the original Hamiltonian. Mathematically, this implies that the following 
expectation value does not vanish: 

(GS|[^^(a;),G]|GS)7^ 0. (235) 
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Note that, in general, we can also substitute [^^, G] for the commutator in Eq. (235), where 
^ is a composite field of and Then, the NG theorem predicts the appearance of 
the corresponding NG mode whose energy vanishes in the long- wavelength limit. 

According to Ref . [22j , the energy of the NG mode obeys a power law of the wave number 
in the long- wavelength limit, and the NG mode is classified as type-I or type-II according 
to whether this power is odd or even, respectively. Nielsen and Chadha formulated the 
following theorem: if the NG mode of type-I is counted once and that of type-II is counted 
twice, then the total number of the NG modes is equal to or greater than the number of the 
symmetry generators that correspond to the spontaneously broken symmetries: 

iVNG = A^i + 2iVii>iVBG, (236) 

where A^i, A^n, Abg are the total number of the type-I NG mode, that of the type-II NG 
mode, that of the symmetry generators whose symmetries are spontaneously broken, respec- 
tively. 

Lorentz invariant theories can have only NG modes with linear dispersion relations. In 
non-relativistic theories, however, there are examples that belong to type-II NG modes. A 
well-known example of type-II NG modes is the Heisenberg ferromagnet, which has a NG 
mode satisfying a quadratic dispersion relation. The authors of Ref. [36] have proved that 
if Gi^ {i = 1, 2, n), constitute a full set of broken charges, and if (GS|[(7i, (7j]]|GS) = for 
any pair (i, j), the relation Ai + An = Abq is satisfied. In Ref. [37], the author suggests that 
if the above commutators are not equal to zero, namely (GS|[(7^, Gj]]|GS) 7^ 0, a type-II NG 
mode appears. 

We apply the above arguments to spinor BECs. In the absence of an external magnetic 
field, the low-energy Hamiltonian in a spinor BEG has the U{1) x 5^0(3) symmetry, rep- 
resenting the global gauge and spin-rotation symmetries. Specifically, in the U{1) 
and SO (3) transformations are given by 

5uii)'i^m{^) = ieifm{^), (237) 

SsO{S)'^m{x) = ^^7w^mK^). U = V z) (238) 

where Q and Q' are infinitesimal parameters. The conserved charges of the t/(l) and 5^0(3) 
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transformations are 

Gu(i) = ^ = j d^¥^{x)^rn{x), (239) 
Gsoiz) ^F' = j d^¥^{x)f^^,i>,^,{x), (240) 

where we drop the infinitesimal parameters. These conserved charge operators are nothing 
but the total number and spin operators. In the presence of an external magnetic field, 
however, the S0{?>) symmetry reduces to 5^0(2), representing the spin-rotation symmetry 
about the direction of the external magnetic field, and the symmetry of the Hamiltonian 
therefore reduces to U{1) x S0{2). Since some or all of these symmetries are spontaneously 
broken in each phase, NG modes are expected to emerge. To find the types of the relevant 
NG modes, it is important to specify the order parameter manifold of each phase. This is 
equivalent to finding the combination of gauge transformation and spin rotations to keep 
the order parameter of each phase unchanged. Such a programme has been carried out in 
Refs. \J\ [211 EH Eni UT] for spin-1 and spin-2 BEG in the absence of an external magnetic 
field, based on the mean-field theory. The number of NG modes and the classification are 
investigated based on the Bogoliubov theory by analyzing the long- wavelength limit of the 
massless modes. 

We first discuss the cases of zero external magnetic field. TABLE I shows the summary 
in spin-1 and spin-2 BEGs. Interestingly, the number of symmetry generators that are 
spontaneously broken A^bg is equal to that of the NG modes N^q in the all phases in 
spin-1 and spin-2 BEGs, namely N^q = Abg, and to the best knowledge of the present 
authors, the equality holds in all systems studied so far. For the spin-1 ferromagnetic BEG, 



it follows from Eqs. (31) and (33) that there exist one type-I and one type-II NG modes, 
which describe the density mode and transverse spin mode, respectively. We note that the 
spin-1 ferromagnetic BEG is similar to the Heisenberg ferromagnet in that the type-II mode 
appears if the following condition holds: 

(GS|[F^F^]|GS) = i(GS|F"|GS) ^0. (241) 

This rule applies also to a spin-2 ferromagnetic BEG. For the other spin-1 and spin-2 phases, 
however, the expectation values of the commutators of the generators are equal to zero. 
Therefore, Ni + An = Abq. In addition, all these phases only have the Bogoliubov modes 
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that are linear and massless. Hence, A^i = A^bg for all phases except for the spin-1 and 
spin-2 ferromagnetic phases. 

The spin-2 nematic phase is special since not all of the massless Bogoliubov modes can be 
interpreted as the NG modes. Both of the uniaxial and biaxial nematic phases have the five 
Bogoliubov modes that are linear and massless, while the number of spontaneously broken 
generators is three for the uniaxial nematic phase and four for the biaxial nematic phase. 
However, since the expectation values of the commutators of the generators are always zero, 
A^i must be equal to A^bg according to Ref. [36]. Therefore, there is one Bogoliubov mode 
that is not the NG mode for the biaxial nematic phase and there are two such modes for 
the uniaxial nematic phases. In each of the uniaxial and biaxial nematic phases, there exists 
the Bogoliubov mode that originates from the fiuctuations with respect to r] and is not 
a NG mode because rj is not related to the symmetry of the Hamiltonian and we cannot 
introduce a conserved current on rj. Here, we note that this Bogoliubov mode can also be 
interpreted as a fiuctuation mode with respect to one of the SO (5) directions other than 
SO (3) directions. This is because we can define the fiuctuation operator around the one of 
SO (5) directions as follows: 

COST] 



V2 



(Sk,2 + - sin7yak,o, (242) 



which is the same as the nematic fiuctuation operator defined in Eq. (126). In addition, one 
of the spin fiuctuation modes is not the NG mode in the uniaxial nematic phase as mentioned 
in Sec. V. B. To show this, we can choose the configuration r] = without loss of generality. 
Then the isotropy groups in the uniaxial spin nematic phase include SO (2)^ which represents 
the rotational symmetry around the z axis. Therefore, the spin mode around the z axis is 
not the NG mode since the symmetry is not broken spontaneously. This is also confirmed 
by the fact that the following expectation value is zero in the uniaxial nematic phase: 

(GS|[l^(x),F^]|GS) = 0, (243) 

where ^^(x) is an arbitrary polynomial in ^^J^(x) and ^^^(x). This implies that the Bogoli- 
ubov mode under consideration is not a NG mode. However, this Bogoliubov mode defined 



in Eq. (137) can be interpreted as the fiuctuation mode around the other direction of S0{5)^ 
since we can define the fiuctuation operator 

^(ak,2-ak-2), (244) 
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TABLE I: Order parameter manifold G/H^ the 


Spin-2 BN 


(C/(l) X 50(2))/Z4 


2 


2 


number of spontaneously broken generators A^bg? 


Spin-2 C 


U{1) X 50(2)7^2 


2 


2 



and that of NG modes A^ng in each phase in the TABLE IL Order parameter manifold GjU, the 
absence of an external magnetic field, where i?4 number of spontaneously broken generators A^bg, 
and T represent the dihedral-four and tetrahedral and that of NG modes A^ng in each phase in the 
groups, respectively. presence of an external magnetic field. 



which is the same as the fluctuation operator defined in Eq. (125). When taken together, 
these massless Bogoliubov modes that are not NG modes can be interpreted as the modes 
related to the S'0(5) symmetry that the mean-field solution in the nematic phase has. This 
is worthy of special mention because it is rather exceptional that a massless mode is not 
a NG mode. The other exceptions are a two-dimensional SU (N) Thirring model [42j and 
a two-dimensional superfluid [35j, both of which have massless modes that are not a NG 
modes because of the Coleman-Hohenberg-Mermin- Wagner theorem [35j . 

Finally, we discuss the cases in the presence of an external magnetic field, which are 
summarized in TABLE Since the symmetry of the Hamiltonian reduces to U{1) x 5^0(2) 
and the commutator of U{1) and S0{2) is zero, we must have A^i + A^n = ^bg- The 
Bogoliubov theory shows that for all phases that we have analyzed, only type-I NG modes 
appear; therefore, A^i = A^bg- To the best knowledge of the present authors, this is always 
the case, that is, only type-I NG modes appear, when (GS|[(7i, Gj]|GS) = for any pair 
(i^j) of broken generators. Hence, it is expected that for arbitrary phases in the presence of 
an external magnetic field, only type-I NG modes emerge whenever spontaneous symmetry 
breaking occurs. 
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VII. SUMMARY AND CONCLUDING REMARKS 



We have derived the Bogohubov spectra and LHY corrections to the GSE, pressure, 
sound velocity, and quantum depletion in the presence of a quadratic Zeeman effect. We 
have shown that in the absence of an external magnetic field, the Bogoliubov eflFective Hamil- 
tonian reduces to the sum of sub-Hamiltonians that can be diagonalized by the standard 
Bogoliubov transformations. A nontrivial example is the nematic phase of a spin-2 BEC 
because we cannot determine the sufficient number of ffuctuation operators that decompose 
the Hamiltonian from the symmetry of the Hamiltonian only. However, because the nematic 
phase features an additional continuous parameter ry, the ffuctuation operator with respect 
to 7] can be constructed and allows the decomposition of the Hamiltonian into the sum of the 
sub-Hamiltonians. Furthermore, because of the r] dependence in the Bogoliubov Hamilto- 
nian, the nematic phase is divided into the uniaxial and biaxial nematic phases by quantum 
ffuctuations even though these phases are degenerate at the mean-ffeld level. Finally, the rj 
dependence removes the degeneracy and causes the phase transition between the uniaxial 
and biaxial nematic phases. 

In the presence of the quadratic Zeeman effect, the magnetic quantum number is no longer 
conserved and different modes are coupled. Therefore, the Hamiltonian must be diagonalized 
by treating the multidimensional Hamiltonian explicitly. By explicitly constructing the 
Bogoliubov transformations, we have obtained the LHY corrections to the GSE, sound 
velocity, and quantum depletion. In the case of a scalar BEC, the LHY corrections depend 
only on the coupling constant and density. For spinor BECs, the LHY corrections also 
depend on the quadratic Zeeman effect except for the ferromagnetic phase. Moreover, taking 
into consideration the fact that the spin-dependent coupling constants are small compared to 
the spin-independent one for the alkali species, the main contribution of the LHY corrections 
arises from the density ffuctuation. The LHY corrections in spinor BECs may be observed 
by making the system strongly correlated by using an optical lattice [43j or by controlling 
an external magnetic ffeld. Since the enhancement of the quantum depletion in an optical 
lattice has already been demonstrated by the MIT group in a scalar ^^Na condensate [16j, 
it should also be possible to apply it to spinor condensates. 

Following the argument by Nielsen and Chadha, we have pointed out the relation between 
the number of symmetry generators that are spontaneously broken A^bg and that of the NG 
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modes A^ng- The NG modes are divided into type-I and type-II according to the dispersion 
laws, and the foUowing inequahty must be satisfied: A^ng = A^i + 2A^n ^ A^bg- In contrast, 
for aU the phases that we have analyzed, it is shown that A^ng = A^bg- The type-II NG 
modes only appear in the spin-1 and spin-2 ferromagnetic phases in the absence of an 
external magnetic field, as in the case of the Heisenberg ferromagnet. Although only type-I 
NG modes emerge for the other phases, nontrivial situations arise for the spin-2 nematic 
phase in the absence of an external magnetic field: there exist the Bogoliubov modes that 
have linear dispersion relations but do not belong to the NG modes. 

We have shown that every configuration of the mean-field ground state that we have 
studied is stable against quantum fiuctuations. However, if the configuration is not a mean- 
field ground state, instabilities set in. For the quadratic spectrum, it implies the Landau 
instability, while for the linear spectrum, it implies the dynamical instability. Furthermore, 
it is expected that the instabilities also occur for the configurations that are not stationary 
solutions of the mean-field theory. We illustrated this for the tetrahedral configuration of the 
cyclic phase, which is not a stationary solution of the mean-field theory for nonzero q and is 
unstable regardless of the sign of q. We hope that our analysis helps the determination of 
the ground-state phase of the spin-2 ^''Rb condensate, which lies in the vicinity of the phase 
boundary between the biaxial nematic and cyclic phases for g < [3]j|. 

States that are unstable in the thermodynamic limit may be stable in a mesoscopic 
regime. For example, for the case of the tetrahedral configuration of the cyclic phase, if the 



lowest wave number is higher than a critical value kc = ^/^M\q\JW ^ the Bogoliubov spectra 
of the spin fiuctuations are positive definite. For the Bogoliubov spectrum of the density 



fiuctuation, the condition of Eq. (228) is needed to be positive. Hence, this configuration can 
be stabilized by the balance between the finite-size and quadratic Zeeman effects. Similar 
arguments can be applied to other configurations. 

Finally, we briefiy comment on several applications related to our analysis. Incorporating 
the trapping effect is an important extension of the present theory. Meanwhile, changing of 
the sign of q with the technique discussed in Ref. [19j also gives rise to nontrivial effect. For 



example, the polar configuration in Eq. (23) becomes unstable and is expected to undergo 



the phase transition to the configuration in Eq. (24) through the dynamical instability when 
the sign of q changes from positive to negative. In this process, the polar direction changes 
from the z axis to a transverse direction, triggering spontaneously symmetry breaking of 
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axisymmetry and dynamical creation of half-quantum vortices around which both the polar 
axis and condensate phase rotate tt [44j. 
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Appendix A: Derivation of the ground-state energies 



1. Polar phase 



In the polar phase, the renormalized and bare coupling constants are related up to the 
second Born approximation by 

_ go + 2g2 _ go + 2g2 go + 2g^^' 1 



+ 



3V ^ 2ek 

k 



c™ + [(c<")= + 2(c™)^]ilj:'l. 



(Al) 



Therefore, the mean-field energy diverges as foUows: 



2 ^ 8M^ 



' 1+2' ^ 



(A2) 



On the other hand, the components of the GSE arising from quantum fluctuations involve 
the k~'^ terms as follows: 
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(A3) 
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Hence, the divergence in Eq. (A2) is canceled with the last terms in Eq. (A3) and the 
convergence of the GSE is ensured. The GSE for c^^^ > in the polar phase is given by 
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The GSE for c^^ < can also be derived similarly. 



2. Broke n-axisymmetry phase 



In the broken-axisymmetry phase, we note the following relations: 



E 



A;2 



.(1) 



+ 



.(1) 



k 



.(1) , ..(1) 



-0 



+ C 



k 



.(1) 



.(l)^2 



M 



v' 1 



Therefore, the mean-field term diverges as 
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On the other hand, in the short- wavelength limit, the components of the GSE arising from 
the quantum fluctuations behave as 



1 



E 



(tk + I - -Bk,/.) + (2£k + "4" - "4" - -Bk,,! - -Ek,()) 



E 



Wk 



,(1) 



+ 



.(1) 



(A8) 



52 



Therefore, the above k ^ terms are canceled with those from the mean-field terms. Thus, 
the convergence of the GSE is ensured and the GSE Eq^ is given by 
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3. Nematic phase 

In the nematic phase, the following relation holds: 
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where we use 



sin^(7/ + 7r/3) + sin^(r/ — 7r/3) + sin^ (r/) = 
sm^(77 + 7r/3) + sm^(77 — tt/S) + sin^ (77) = 



Therefore, the mean-field energy is rewritten as foUows: 
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On the other hand, the components of the GSE arising from the quantum fluctuations 
diverge as foUows: 

+ nci'^ + ncf - E^^a) + (^k + ncf\r] + ^/3) - E^j^j 
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As expected, these divergences are canceled out with the last three terms of Eq. (A13), and 
the GSE in the nematic phase Eq is given by 
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We note that the above renormalization procedure can be used in the cases of the uniaxial 
and biaxial nematic phases. To see this, we focus on the uniaxial spin nematic phase. As 
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r] = in the uniaxial nematic phase, the relation (A13) is rewritten as 
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The mean-field term involves the k ^ terms as follows: 
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On the other hand, the components of the GSE arising from the quantum fluctuations 
diverge as foUows: 
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As in the case of g = 0, these divergences are canceled out with the last three terms of Eq. 



(A17) and the finite GSE is obtained. The same holds in the biaxial nematic phase. 



4. Cyclic phase 



In the cyclic phase, we consider the following relations: 
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Hence, the mean-field terms behave as 
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At short- wavelength hmit, the components of the GSE stemming from the quantum fluctu- 
ations diverge as follows: 
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These divergences are canceled out with fc-dependent terms of Eq. (A21). The GSE 
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Appendix B: Lists of equation numbers and symbols 



In this appendix, we list the properties in each phase of spin-1 and spin-2 BECs in Table 



III, equation numbers of various physical quantities in Table IV, and symbols in Table M 



TABLE III: Possible ground-state phases, order parameters (^^, magnetization, and spin-singlet 
pair amplitude defined by (s-) = | X]m(~^)^CmC-m of spin-1 and spin-2 BECs. Here, F, P, BA, 
N, BN, UN, and C stand for ferromagnetic, polar, broken- axisymmetry, nematic, biaxial nematic, 
uniaxial nematic, and cyclic phases, respectively In the spin-2 nematic phase, 77 is an additional 
continuous parameter that represents the degeneracy of the uniaxial and biaxial nematic phases. 
In the spin-2 broken axisymmetry phase, + (— ) sign corresponds to the case of < (> 0). 



Phase 


order parameter 


{{) 


2 X 


Spin-1 F 




(1,0,0) 




in = 1 





P 




(0,1,0) 




(f) = 


1 


P' 




(l/AO,l/v^) 




(f) = o 


1 


BA 


(\/l/4 


+ <z/(8ncW), ^ _ q/(4ncW), ^1/4 + g/(8nc«)) 


(r) = 


^Jl - (g/2ncW)2 


-g/(2n4^)) 


Spin-2 F 




(1,0,0,0,0) 




in = 2 





N 




(sin77/\/2, 0, cos 77, 0, smri/\/2) 




(f) = 


1 


BN 




(l/AO,0,0,l/v^) 




(f) = 


1 


UN 




(0,0,1,0,0) 




(f) = 


1 


C 


(\/l/4- 


q/(2ncf ), 0, ^J 1/2 + <//(ncf ), 0, ^ 1/4 - q/(2ncf )) 




(f) = 


2g/(ncf) 


M 


( 


^ 1/3 - q/(3ncf ^), 0, 0, ^ 2/3 + q/{3ncf^), o) 


in- 


= -4g/(3n4^^) 







or 


(O, ^J2/3 + g/(3ncp)),0,0, y^l/3 - q/i^ncf^)) 


in 


= 4g/(3ncf^) 





BA 




(iba, b, c, , 6, ±a) 


{n = 


4b{^y3/2c±a) 


2(a2-62) + c2 
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TABLE IV: Equation numbers of various physical quantities. 



Phase 


Bogoliubov spectra 


GSE 


Pressure 


Sound velocity 


Depletion 


Spin-1 F 

P 

P' 

BA 


( 
( 


47 
78 


( 

), ( 
), ( 


33 
48 
67 
86 


) 
) 
) 
) 


( 
( 
( 
( 


36 
50 
68 
89 


) 
) 
) 
) 


( 
( 
( 
( 


37 
52 
69 
92 


) 
) 
) 
) 


( 
( 
( 
( 


38 
53 
70 
93 


) 
) 
) 
) 


( 
( 
( 
( 


39) 
55) 
71) 
95) 


Spin-2 F 
N 

BN 
UN 
C 


( 


184 


( 
( 
( 

), ( 


134 
156 
170 
188) 


( 

)-( 
)-( 
)-( 

, ( 


115 
138 
159 
172 
195 


) 
) 
) 
) 
) 


( 
( 
( 
( 
( 


118 
141 
160 
173 
198 


) 
) 
) 
) 
) 


( 
( 
( 
( 
( 


119 
142 
161 
174 
202 


) 
) 
) 
) 
) 


( 
( 
( 
( 
( 


120 
143 
162 
175 
203 


) 
) 
) 
) 
) 


( 
( 
( 
( 
( 


121) 
144) 
163) 
176) 
206) 
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TABLE V: List of symbols. Note that has the dimension of the couphng constant and U is 
dimensionless. 



Common 


q {gfiBB)yEM 


Spin-1 BEC 


4 (90 + 252)/3 
c? (92 - 50)73 
gV(4n2cW) 

h 97rfl)-i 
*2 -57(«4'^) 

*3 (?V(n4'^+n4'))2 


Spin-2 BEC 


4') (452+354)77 

cf^ (54-52)77 

4'^ (750-1052+354)735 

4^^(??) 44^^sin2r7-4^^ 
4'^ 4gV (71242)) 

I //|(z) -i 

U -2,q/{n\c\ ' -c^'y) - I 
h -457(n|4')|)-l 
te ^7(n|34')-4')|)-l 
t, q/in\cf^\)-l 
ts \q\/(2ncf^) 
t, (2q)V{nc^^^r 



[1] N. N. Bogoliubov, J. Phys. (U.S.S.R.) 11, 23 (1947). 

[2] T. D. Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957). 

[3] T. D. Lee, K. Huang, and C. N. Yang, Phys. Rev. 106, 1135 (1957). 

[4] T. T. Wu, Phys. Rev. 115, 1390 (1959). 

[5] N. M. Hugenholtz and D. Pines, Phys. Rev. 116, 489 (1959). 

[6] K. Sawada, Phys. Rev. 116, 1344 (1959). 



60 



[7] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[8] T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 (1998). 

[9] W. -J. Huang and S. -C. Gou, Phys. Rev. A 59, 4608 (1999). 
[10] M. Ueda, Phys. Rev. A 63, 013601 (2000). 
[11] P. Szepfalusy and G. Szirmai, Phys. Rev. A 65, 043602 (2002). 
[12] K. Murata, H. Saito, and M. Ueda, Phys. Rev. A 75, 013607 (2007). 
[13] J. Ruostekoski and Z. Button, Phys. Rev. A 76, 063607 (2007). 
[14] J.-P. Martikainen and K.-A. Suominen, J. Phys. B 34, 4091 (2001). 
[15] M. Ueda and M. Koashi, Phys. Rev. A 65, 063602 (2002). 

[16] K. Xu, Y. Liu, D. E. Miller, J. K. Chin, W.Setiawan, and W. Ketterle, Phys. Rev. Lett. 96, 
180405 (2006). 

[17] S. B. Papp, J. M. Pino, R. J. Wild, S. Ronen, C. E. Wieman, D. S. Jin, and E. A. Cornell, 
Phys. Rev. Lett. 101, 135301 (2008). 

[18] A. Altmeyer, S. Riedl, C. Kohstall, M. J. Wright, R. Geursen, M. Bartenstein, C. Chin, J. H. 
Denschlag, R. Grimm, Phys. Rev. Lett. 98, 040401 (2007). 

[19] S. R. Leslie, J. Guzman, M. Vengalattore, J. D. Sau, M. L. Cohen, and D. M. Stamper-Kurn, 
Phys. Rev. A 79, 043631 (2009). 

[20] A. M. Turner, R. Barnett, E. Dernier, and A. Vishwanath, Phys. Rev. Lett. 98, 190404 (2007). 

[21] J. L. Song, G. W. Semenoff, and F. Zhou, Phys. Rev. Lett. 98, 160408 (2007). 

[22] H. B. Nielsen and S. Chadha, Nucl. Phys. B 105, 445 (1976). 

[23] M. Koashi and M. Ueda, Phys. Rev. Lett. 84, 1066 (2000). 

[24] C.V. Ciobanu, S.-K. Yip, and T.-L. Ho, Phys. Rev. A 61, 033607 (2000). 

[25] R. Barnett, A. Turner, and E. Dernier, Phys. Rev. Lett. 97, 180412 (2006). 

[26] See, e.g., C. J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases, (Cam- 
bridge University Press, Cambridge, 2008). 

[27] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Miesner, A. P. Chikkatur, and W. Ketterle, 
Nature (London) 396, 345 (1998). 

[28] H. Saito and M. Ueda, Phys. Rev. A 72, 053628 (2005). 

[29] N. D. Mermin, Phys. Rev. A 9, 868 (1974). 

[30] A. Widera, F. Gerbier, S. Foiling, T. Gericke, O. Mandel, and L Bloch, New J. Phys. 8 152 
(2006). 

61 



[31] S. Tojo, T. Hayashi, T. Tanabe, T. Hirano, Y. Kawaguchi, H. Saito, and M. Ueda, Phys. Rev. 
A 80, 042704 (2009). 

[32] Since, the spin fluctuation around the z axis vanishes when r/ = 0, we cannot identify Eq. 



(125) as the spin fluctuation. However, we can interpret Eq. (125) as the fluctuation operator 

with respect to one of the 50(5) directions as discussed in Sec. VI. 
[33] S. Uchino, T. Otsuka, and M. Ueda, Phys. Rev. A 78, 023609 (2008). 
[34] T. Brauner, Symmetry 2, 609 (2010). 

[35] N. Nagaosa, Quantum Field Theory in Condensed Matter Physics (Springer- Verlag, Berhn, 
1999). 

[36] T. Schafer, D. T. Son, M. A. Stephanov, D. Toublan, J. J. M. Verbaarschot, Phys. Lett. B 

522, 67 (2001). 
[37] T. Brauner, Phys. Rev. D 72, 076002 (2005). 
[38] F. Zhou, Phys. Rev. Lett. 87, 080401 (2001). 

[39] H. Makela, Y. Zhang, and K.-A. Suominen, J. Phys. A: Math. Gen. 36, 8555 (2003). 

[40] M. Kobayashi, Y. Kawaguchi, M. Nitta, and M. Ueda, Phys. Rev. Lett. 103, 115301 (2009). 

[41] G. W. Semenoff and F. Zhou, Phys. Rev. Lett. 98, 100401 (2007). 

[42] E. Witten, Nucl. Phys. B 145, 110 (1978). 

[43] J. L. Song and F. Zhou, Phys. Rev. A 77 033628 (2008). 

[44] S. Uchino, M. Kobayashi, and M. Ueda (to be pubhshed). 



62 



